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ON  THE  KERNEL  FUNCTION  OF  THE  INTEGRAL  EQUATION  RELATING  THE  LIFT  AND 
DOWNWASH  DISTRIBUTIONS  OF  OSCILLATING  FINITE  WINGS 

IN  SUBSONIC  FLOW  ‘ 

By  Chables  E.  Watkins,  Habbv  L.  Runyan,  and  Donald  S.  Woolston 


SUMMARY 

This  report  treats  the  kernel  function  of  an  integral  equation 
UuU  relates  a  known  or  prescribed  doumwash  distribution  to  an 
unknown  lift  distribution  for  a  harmonically  oscillating  finite 
wing  in  compressible  subsonic  flow.  The  kernel  function  is 
reduced  to  a  form  that  can  be  accurately  evaluated  by  separating 
the  kernel  function  into  two  parts:  a  part  in  which  the  singular¬ 
ities  are  isolated  and  analytically  expressed  and  a  nonsingular 
part  which  may  be  tabulated.  The  form  of  the  kernel  function 
for  the  sonic  case  (Mach  number  of  1)  is  treated  separately.  In 
addition,  results  for  the  special  cases  of  Mach  number  of  0 
(incompressible  case)  and  frequency  of  0  (steady  case)  are  given. 

The  derivation  of  the  integral  equation  which  involves  this 
kernel  function,  originally  performed  elsewhere  (see,  for  example, 
NACA  Technical  Memorandum  979),  is  reproduced  as  an 
appendix.  Another  appendix  gives  the  reduction  of  the  form  of 
the  kernel  function  obtained  herein  for  the  three-dimensional 
case  to  a  known  result  of  Possio  for  two-dimensional  flow.  A 
third  appendix  contains  some  remarks  on  the  evaluation  of  the 
kernel  function,  and  a  fourth  appendix  presents  an  alternate 
form  of  expression  for  the  kernel  function. 

INTRODUCTION 

The  analytical  determination  of  air  forces  on  oscillating 
wings  in  subsonic  flow  has  been  a  continuing  problem  for  the 
past  30  years.  Throughout  the  flrst  and  greater  part  of 
this  time,  efforts  were  directed  mainly  toward  the  determina¬ 
tion  of  forces  on  wings  in  incompressible  flow.  These  efforts 
have  led  to  important  closed-form  solutions  for  rigid  wings 
in  two-dimensional  flow  (ref.  1),  to  solutions  in  terms  of 
series  of  Legendre  functions  for  distorting  wings  of  circular 
plan  form  (refs.  2  and  3),  and  to  many  approximate,  yet 
useful,  results  for  wings  of  elliptic,  rectangular,  and  tri¬ 
angular  plan  form  (see,  for  example,  refs.  4  to  12). 

Although  these  results  for  incompressible  flow  play  a 
highly  significant  role  in  applications  of  unsteady  aerody¬ 
namic  theory,  the  advent  of  higher  and  higher  speed  aircraft 
during  the  last  15  years  has  brought  a  growing  need  for 
knowledge  of  the  effect  that  the  compressibility  of  air  might 
have  on  unsteady  air  forces,  or  for  analytically  derived  un¬ 
steady  air  forces  based  on  a  compressible  medium.  The 
transition  to  results  for  a  compressible  fluid  from  those  for 
an  incompressible  fluid  is  not  likely  to  be  accomplished  by 
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applications  of  sunplc  transformations  or  correction  factors, 
such  as  tlic  well-known  Prandtl-GIauert  factor  for  steady 
flow.  This  difficulty  is  associatcrl  with  the  fact  that  the  time 
requiretl  for  signals  arising  at  one  point  in  the  medium  to 
reach  other  points  gives  rise  not  only  to  changes  in  magni¬ 
tudes  of  forces  but  also  to  additional  phase  lags  between 
instantaneous  positions,  velocities,  and  accelerations  of  the 
wing  and  tlic  corresponding  instantaneous  forces  associated 
with  these  quantities.  In  order  to  obtain  results  for  the 
compressible  case,  it  therefore  appears  necessary  to  deal 
directly  with  the  boundary-value  problem  for  this  case. 

The  boundary-value  problem  for  a  two-dimensional  wing 
in  compressible  flow  lias  been  successfully  attacked  from  two 
points  of  view.  First,  by  consideration  of  an  acceleration  or 
pressure  potential,  Possio  (ref.  13)  reduced  the  problem  to  that 
of  an  integral  cqtialion  relating  a  prescribed  downwaah  dis¬ 
tribution  to  an  unknown  lift  distribution.  The  kernel  of  this 
integral  equation,  which  is  a  rather  abstruse  function,  was 
reduced  to  a  form  that,  c.\cept  at  singular  points,  could  be 
evaluated.  Schwarz  (ref.  14)  later  isolated  and  determined 
the  analytic  behavior  of  the  singular  points  of  Possio’s  results 
and  made  fairly  extensive  tables  of  the  kernel  function. 
These  tabular  values  were  use<l  by  various  investigators 
(for  examples,  refs.  15  and  16)  to  obtain,  by  numerical 
procedures,  initial  tables  of  force  and  moment  coefficients 
for  oscillating  wings  in  compressible  subsonic  flow. 

The  second  successful  approach  to  the  solution  of  the 
boundary-value  problem  for  a  two-dimensional  wing  (see 
refs.  17  to  19)  is  achieved  by  a  transformation  to  elliptic 
coordinates  followed  by  a  separation  of  variables  that  reduces 
the  boundary-value  problem  from  one  in  partial-differential 
equations  to  one  in  ordinary  differential  equations  of  the 
Mathieu  typo.  The  solutions  turn  out  os  infinite  series  in 
terms  of  ^Iathicu  functions.  Numerical  results  obtained 
recently  by  this  procedure  agree  with  results  previously  ob¬ 
tained  by  the  numerical  procedures  using  the  kernel  func¬ 
tion  (see,  for  example,  ref.  20). 

With  regard  to  boundary-value  problems  for  finite  nhigs 
in  compressible  flow,  it  appears  that  the  procedure  of  sepa¬ 
ration  of  variables  could  be  a  feasible  approach  only  for 
wings  of  very  special  plan  forms  such  as  a  circle  or  an  ellipse. 
In  any  case,  the  development  of  the  appropriate  mathe¬ 
matical  functions  for  a  particular  plan  form  would  become 
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highly  involvotl.  On  the  other  hand,  it  appears  that  approxi¬ 
mate  procedures  similar  to  those  used  for  two-dimensional 
wings  might  afford  an  approach  to  solutions  of  these  prob¬ 
lems  which,  though  laborious,  might  be  handled  by  routine 
numerical  methods. 

The  kOTiel  function  of  the  integral  equation  relating  pres¬ 
sure  and  downwash  for  the  three-dimensional  case  appears 
as  an  improper  integral.  The  purpose  of  this  report  is  to 
treat  and  discuss  this  kernel  function.  The  improper  integral 
is  redui^  to  a  form  that  can  be  accurately  evaluated  by 
numerical  procedures.  The  form  and  order  of  all  its  singular¬ 
ities  are  determined  and  an  expression  for  the  kernel  function 
is  derived  in  which  the  singularities  are  isolated.  Special 
forms  of  the  kernel  for  the  sonic  case  {M=  1) ,  the  incompres¬ 
sible  case  (M=0),  and  the  steady  case  (ii;=0)  are  presented. 
A  series  expansion  in  powers  of  the  reduced-frequency  param¬ 
eter  k  is  developed. 

The  availability  of  the  kernel  in  a  form  which  can  be 
rapidly  evaluated-  makes  possible  the  use  of  numerical  pro¬ 
cedures,  similar  to  those  used  in  the  two-dimensional  case, 
to  obtain  aerodynamic  forces  for  finite  vdngs. 

SYMBOLS 


c 

velocity  of  sound 

Hankel  functions  of  second  kind  of  zero 

and  first  order,  respectively 
modified  Bessel  functions  of  first  kind  of 
zero  and  first  order,  respectively 

Bessel  function  of  first  kind  of  zero  order 
modified  Bessel  functions  of  second  kind  of 

J* 

zero  and  first  order,  respectively 

Kixt,y9) 

kernel  ftmetion  of  integral  equation 

E'(xo,yit) 

singular  part  of  K(xb,yo) 

k 

reduced-frequency  parameter,  la/V 

modified  Struve  functions  of  zero  and  first 

order,  respectively 

unknown  lift  distribution 

1 

reference  length 

M 

Mach  number,  V/e 

P 

pressure 

S 

region  of  zy-plane. occupied  by  wing 

t 

time 

V 

forward  velocity  of  wing 

w{x,y) 

amplitude  fimction  of  prescribed  downwash. 

v>(x,y,t)=e'^(x,y) 

Cartesian  coordinates 

ih=y-v 

y 

Euler’s  constant 

«=Vxo’+/S*yo’ 

* 

velocity  potential 

* 

acceleration  potential 

P 

fluid  density 

a 

circular  frequency  of  oscillation 

'S=ulV^ 


ANALYSIS 

INTBGRAL  EQUATION  AND  OMGINAL  FORM  OF  RERNRL  FUNCTION 

The  main  purpose  of  this  analysis  is  to  treat  the  kernel 
function  of  an  integral  equation  that  relates  a  known  or 
prescribed  downwash  distribution  to  an  unknown  lift  dis¬ 
tribution  for  a  harmonically  oscillating  finite  wing  in  com¬ 
pressible  subsonic  flow.  The  integral  equation  referred  to 
can  be  obtained  by  employing  the  Prandtl  acceleration 
potential  to  treat  linearized  boundary-value  problems  for 
oscillating  finite  wings  by  means  of  doublet  distributions. 
Derivation  of  this  integral  equation  from  the  linearized 
boundary-value  problem  for  a  wing  is  a  preliminary  task 
that  has  been  done  elsewhere  (see,  for  example,  ref.  21),  but 
it  is  reproduced  herein  as  an  appendix  for  the  sake  of  com¬ 
pleteness. 

In  keeping  with  the  concepts  of  linear  theory,  the  wing  is 
considered  a  plane  impenetrable  surface  S  which  lies  neaHy 
in  the  zy-plane  as  indicated  in  sketch  1 : 


/ 


The  x,y,z  coordinate  system  and  the  surface  S  are  assumed  to 
move  in  the  negative  z-direction  at  a  uniform  velocity  V. 

In  terms  of  these  coordinates,  the  integral  equation  may  be 
formally  written  as 

where  ^(z,y)  is  the  amplitude  function  of  the  prescribed 
downwash,  it(a%,yo)=K(i— {,  y—v)  ia  the  kernel  function 
and  physically  represents  the  contribution  to  downwash  at 
a  field  point  {x,y)  due  to  a  pulsating  pressure  doublet  of  unit 
strength  located  at  any  point  and  L  ((,v)  is  the  unknown 

lift  distribution  or  local  doublet  strength. 
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The  kernel  function  may  be  mathematically  defined  by  the 
following  improper  integral  eglpression  (see  eq.  (A12), 
appendix  A) : 


VX*+^o'+/3*2* 


dX  (2) 


where  M  is  Mach  number,  /3=Vl— M*,  u=u/V^,  u  is  the 
circular  frequency  of  oscillation,  V  is  the  velocity,  and 
X  is  the  variable  of  integration.  Evaluation  of  this  integral 
constitutes  a  main  difficulty  in  obtaining  aerodynamic 
coefficients  for  oscillating  finite  wings  in  compressible  flow. 
The  present  analysis  is  therefore  devoted  to  reducing  it  to  a 
form  that  can  be  accurately  evaluated  by  numerical  pro* 
cedures  combined  with  the  use  of  tables  of  certain  tabulated 
functions.  The  form  and  order  of  all  its  angularities  are 
detennined,  and  an  expression  for  the  kernel  function  is 
derived  in  which  the  sincerities  are  isolated. 

eaDUcnoN  or  n»  eunbl  roNcnoN 

In  considering  the  reduction  of  the  kernel  function 
K(xo,yo),  the  int^ral  involved  can,  for  convenience,  bo 
written  as  the  sum  of  two  integrals,  namely 


•t,  «c(x-vVS3+^ 

V5?+? 


X 


g-i5(x+vVx>+f») 

•tt 

V5?+? 


dX+ 


dX 


(3) 


Therefore, 


a*  -! 


1^-lim  ^e~y  (F)=:lim  ^ <  >'  (F.+f,)  (4) 


where 


and 


H. 

ri,  «C(x-MV^>+r>) 


(5) 


(6) 


and  where  r= 

The  integrals  Ft  and  Ft  are  treated  separately  in  succeeding 
sections.  The  final  forms  are  given  in  equations  (15)  and 
(19),  respectively. 

Ivalnation  of  Fi. — ^The  integral  Fi  can  be  converted  to  a 
form  that  can  be  more  easily  handled  by  writing 


®  — c-r  -r~ 


and  introducing  the  following  relation  (see  p.  416  of  ref.  22) 
«-BMVxi+S  p- ^  , 


■X 

jMim 


ij^  j,(rx) 


VM*S*-2** 


TdT 


(7) 


In  the  first  integral  of  these  last  two  intograb,  make  tlie 
substitution 

Vr^Af»S*=T 

and  in  the  second  integral  make  tlie  substitution 


Then 


Jo(XvV^+3J^  dr- 
/•MS  ^ _ 

i  I  dr 


(8) 


dx  [J]“  J,(xV?7^l^)  dr- 

dr] 


(11) 


or 


(It  b  of  interest  to  note,  in  the  e.xpre8siou  on  the  left  of  eq. 
(8),  that  X  and  r  appear  in  the  same  manner.  The  roles  of 
thexs  two  quantities  could,  therefore,  be  interchanged  in 
the  expression  on  the  right.) 

With  use  of  equation  (8),  the  equation  for  Fi  can  be 
written  as 


(8) 


ay 

0/. 


Changing  the  order  of  integration  in  each  integral  (which  b 
a  legitimate  step  because  the  int^rahds  inCved  satisfy 
the  continuity  conditions  required  for  such  operatioiis)  leads 
to  the  following  expression  for  Ft: 

F,=’f‘  e-"  dr  «-<»  Jo  (xV?+Af*^)  dxj- 

The  integrab  within  the  brackets  in  equation  (10)  may  be 
evaluated  from  tables  of  Fourier  or  Laplace  transforms  as 
(see,  for  example,  pair  no.  55  of  appendix  III  of  ref.  23) 

Jo  (X^M*5»)  d\=^y^ 

so  that 

„  f  e-"  j  f"®  j 

Jo  V?+^ 

The  first  integral  in  equation  (11)  can  be  written  as 


4/ 


J 


f’-J^Z=dr=f“e-»‘^^dff-^r'‘e-^‘"*df  (11a) 

Jo  Vt*— S’3*  Jo  2J_r/» 
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The  first  integral  on  the  right  of  equation  (11a)  is  given  on 
page  181  of  reference  22  as 

J  g-«sr  c<»h » de=Ko(0Zr) 

where  is  the  modified  Bessel  function  of  the  second  kind 
of  zero  order.  The  second  integral  on  the  right  of  equation 
(11a)  is  given  on  page  338  of  reference  22  as 

'  <^=-Y  [/o(l»i3r)-7V/*Sr)J 

where  !o  is  the  modified  Bessel  function  of  the  first  kind  of 
zero  order  and  La  is  the  modified  Struve  function  of  zero 
order.  Then,  the  first  integral  of  equation  (11)  can  bo 
written  as 

Jo  (fT=go(/3wr)— ^  [/a(j3ar)—£a03ur)] 

(12) 

Note  that  the  end  result  indicated  in  equation  (12)  is  in* 
dependent  of  Mach  number.  The  second  integral  in  equa¬ 
tion  (11)  may  be  written  in  another  form  as 

TMS  g-irr 

I  ,-jv  .j;  dT=  I  - ?=f=i - (13) 

Jo  J®  VH*r* 

This  integral  has  not  been  reduced  to  closed  form;  however, 
it  is  nonsingular  and  can  be  readily  handled  by  numerical 
methods. _ 

(??)  ®W(ainh9-M  cosh  «)«*»•*( 


Combining  equations  (12)  and  (13)  gives  the  following 
c.xprcssion  for  Fi : 

f i=ii:o  Vy7+?)-f  [lo  (9  Vi^)- 


i.(pvs?+?)]-/; 


ViT? 


By  performing  the  dilferentiatious  indicated  in  equation  (4), 
there  is  obtained  for  the  first  part  of  equation  (4)  the  follow¬ 
ing  expression: 

j  J|y.l  dr}  (15) 

All  terms  of  this  e.\pression  other  than  the  integral  may 
be  evaluated  at  small  intervals  of  ya  from  existing  tables, 
except  at  yo=0  where  the  function  is  singular.  The  integral 
is  well  behaved  and  can  be  accurately  evaluated  by  numerical 
or  approximate  procedures.  The  t3rpe  and  order  of  the 
singularities  at  yo=0  sre  discussed  in  a  later  section. 

Evaluation  of  Fj. — In  order  to  reduce  the  integral  Fa, 
equation  (6),  it  is  convenient  to  make  the  substitution 


X=!r  siinh  $ 


so  that 


J’slnh~‘5 


f  fir  (dab  «-Mooib  «)  ( 


Noting  that  z  appears  only  in  r  and  performing  the  differen¬ 
tiations  indicated  in  equation  (4)  yidds 

dg - ^  eC(i.-Arv»+w) 

iA)VV+fi*yo’ 


[fi*  cosh  9— (cosh  $—M  sinh  9))  eO»l»«l(diih»-Moodi»>  de 


-  ff  1  J  —e-BMstn)]  — r’*”"  cosh  9 e®*''*' «i9 }  (jg) 

\  /JVoVa.*+W  ^  J 


or,  by  reverting  completely  to  Cartesian  coordinates  through 
equation  (16),  there  is  obtained 

V^A-o  ^  yaSxa*+(h^ 

(19) 

This  expression  vanishes,  as  it  should,  for  Zo=0  and,  like 
that  in  equation  (15),  has  singularities  at  yo=0  which,  also, 
will  be  handled  in  a  later  section.  The  int^;ral  that  remains, 


like  the  integral  remaining  in  equation  (15),  is  nonsingular 
and  simple  in  form  and  can  be  readily  evaluated  by  numerical 
procedures. 

Expression  for  the  kernel  in  terms  of  nondimensional 
length  variables. — Equations  (15)  and  (19)  can  now  be 
combined  to  give  a  reduced  form  of  the  kernel  function 
X(xa,yt).  However,  in  application,  the  variables  xo  and  ya 
are  employed,  for  convenience,  in  nondimensional  form. 
This  is  accomplished  by  considering  these  variables  in  a 
new  sense  to  mean  that  they  have  been  referred  to  some 
chosen  length  I  and  by  introducing  the  reduced-frequency 
parameter  k=lulV.  The  variables  will  be  used  in  this  new 
sense  throughout  the  remainder  of  the  report.  The  kernel 
can  be  written  in  terms  of  these  nondimensional  variables  as 
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M(ky^\ 


f  ' 

An  alternate  and  perhaps  more  desirable  form  of  expression 
for  the  kernel  function  is  given’ in  appendix  D. 

Note  that  this  expression  for  K(Xi,yt)  can  be  considered  as 
a  function  of  only  three  parameters,  namely,  klyol,  Icx^,  end  M. 
To  be  more  specie,  the  first  two  terms  are  functions  only  of 
k|yol ;  the  next  two  terms  are  functions  of  il:|yo|  and  M;  and 
the  last  two  terms  are  fimctiona  of  klyoli  end  M. 

Equation  (20)  constitutes  the  principal  result  of  this  report. 
Some  partial  checks  as  to  its  correctness  are:  (1)  For  k=0, 
it  reduces,  as  discussed  subsequently,  to  the  downwash  of  a 
pressure  doublet  in  steady  fiow  and  (2)  an  integration  with 
regard  to  the  y-direction  between  the  limits  to  + «» 
yields  Possio’s  result  for  the  two-dimensional  case.  This 
integration  is  carried  out  in  appendix  B.  Other  special  forms 
of  the  kernel  function  for  M=  1,  M=0,  and  &=0  are  derived 
in  subsequent  sections.  A  power  series  expansion  of  the 
kernel  which  is  applicable  for  certain  ranges  of  tbe  parameters 
f’lVtit  hxt,  and  M  is  presented.  In  the  section  immediately 
following,  the  orders  and  of  the  singularities  of  the 
kem^  function  are  discussed. 

nncrasioN  or  tbs  siNOULAarms  or  the  aaaNSL  ruNcnoN 

As  previously  indicated,  the  kernel  function  becomes  singu¬ 
lar  or  indetenninate  at  yo=0.  The  forms  that  the  kernel 
function  takes  when  it  becomes  singular  are  of  particular 
importance  in  applications  to  lifting  surface  theory.  It  is 
therefore  desirable  to  extract  and  treai£)  the  sin^arities 
separately. 

This  extraction  can  be  conveniently  made  by  considering 
tbe  value  of  X(2^,ya),  equation  (20),  at  points  on  tbe  semi¬ 
circumference  of  a  small  ellipse  (see  sketch  2),  the  polar 
equation  of  which  may  be  written  as 


2^=se  sin  9 

Vt^coae  j 


where,  because  of  the  symmetry  of  K(xo,ya)  with  respect  to 
ye,  only  tbe  limits  —  ■■/2Si0$r/2  need  be  examined.  Note 
that  in  these  equations  values  of  9  in  the  range  —  r/2  $0<O 
correspond  to  field  points  ahead  of  or  upstream  from  the 
doublet  position  and  values  of  9  in  the  range  0<9:St/2,  to 
field  points  behind  or  downstream  from  the  doublet  position. 
In  particular,  8=xl2  corresponds  to  points  directly  behind  or 
in  the  wake  of  the  doublet. 


Sketch  2. 


After  substituting  these  expressions  for  Xt  and  into 
equation  (20),  the  results  may  be  written  as 

trt  iis—lS**"'***'**  r  !:«  cos  9  ,, /if  cos  9\ 

\ - T~ 

[/,  (^*)> 


[/,  (^*)]+ 

<Uf>  «M  •  !*•  i»ia  f-if) 

«  s*  ,  tit  cos  9  1 

-na — —M — +— flT- «  '  ^ 


P  -r*  ^ 


ik  , 

irj.  * 


?V) 

dxj-  (22) 


With  the  use  of  the  follovring  series  expressions  for  £i(r)  and 
[/i(s)  —Lx{z)]  (which  can  be  obtained  from  ref.  22 — for  Kx, 
see  p.  80;  for  see  p.  77;  and  for  £],  see  p.  320) : 

where  r  is  Euler’s  constant  (t3=0.5772157),  and 

•  •  •  S4) 

it  is  found  that  for  vanishingly  small  values  of  f  the  limiting 
value  of  the  expression  for  K(€,8)  in  equation  (22)  is  for 
M<1 

f  '  ^  ii(l— sin  9)  , 

*(*’®)“^r“\?(l-8iii9)+T“2*o«  2(1-S)  “ 
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where  0(i")  represents  terms  of  order  c*  for  n  ^  1 .  Expressed 
in  terms  of  lo  and  j/o.  equation  (25)  becomes 


K(X(„yo)‘ 


‘-iri 


—(x<,+\i(i‘+0*yo^  .  ik 
J/oVio*+^o’ 


Examination  of  equation  (25)  shows  that  the  kernel  function 
K(*,9)  has  singularities  with  respect  to  *=-)fx7+Wy?  as 
follows: 

T*  -^UM+logt]  (27) 


where,  from  equation  (25), 


i8*(l+fr  fl); 
cos’  8 


....  ,  iKl— aintf) 
M9)=^og 


k  cos’tf 


2(l-A0(l+sin«) 


Although  of  no  particular  significance  in  applications,  it  is 
of  interest  to  note  that  the  quantities  Ji  and  /i  each  have 

minimum  values  and  |/jUi«=log 

fs:— «r/2,  which  corresponds  to  points  directly  ahead  of  the 
doublet  position;  and,  as  9  increases  from  —  to  +v/^,  the 
values  of  these  quantities  continuously  increase  from  these 
minimum  values  to  infinite  quantities  as  follows: 

/.  ( =lim  =li'n  V 

ico8*(j-«)  ►  (29) 

/.  (I )  =lim  log - r  -- 

V2/  2(l-iW)[l+8in(|-«)] 

..  ,  *«’  I 

•>lim  log  575 - jra 

^4(1— Af)|  J 

Thus  K(xg,yo)  is  singular  for  9=t/2  even  when  the  distance 
«  from  the  doublet  is  not  necessarily  of  zero  order.  This 
implies  that  the  doublet  produces  a  wake  of  discontinuous 
downwash  that  extends  downstream  from  the  doublet 
position  to  infinity. 

With  knowledge  of  the  singularities  involved  in  the  kernel 
function  K(xo,yo),  an  expression  can  be  written  in  which  the 
kernel  is  separated  into  a  singular  part  and  a  nonsingular 
part  (as  was  done  by  Schwarz,  ref.  14,  for  the  two-dimen¬ 
sional  case)  as  foUows 

K(xt,yo)  ■  [K(xt,yo)  -  K'  (r6,yo)l + if'(xb,J/o)  (30) 

where  fiTCzbiVo)  is  defined  in  equation  (20)  or  (22)  and 


K'(t  .  ik 

°  L  yoV*o’+/3V 

^  Xt—M\Xo^+&^yB^  P  I  Jo)1 

2^’  2  2(1 -A/)  J 


2^’  y'iTWW  2  2(1 -A/)  J 

or  in  terms  of  <  and  9,  introduced  by  equations  (21), 

<»> 

The  term  lK(xB,yB)—K'(xB,ya)]  in  equation  (30)  is  a  continuous 
function  for  all  values  of  k,  Zo,  and  Vo  and  for  values  of  Af  in 
the  range  of  0  ^  Af  ^  1.  The  term  K'(xo,ya)  is  discontinuous 
at  the  doublet  position  (zo=0,  Vo=0)  and  at  all  points  in 
the  wake  (•ro>0,yo=0).  It  is  to  be  noted,  however,  that 
each  term  of  Ar'(zo,yo)  possesses  a  simple  indefinite  integral 
with  respect  to  y#  or  with  respect  to  i|=y— Voi  s  fset  that 
may  be  useful  in  some  numerical  applications.  The  manner 
in  which  these  integrals  are  to  be  evaluated  is  indicated  in 
a  subsequent  section  that  deals  with  steady  flow.  The 
limiting  values  at  ys=0  of  lK(xB,ya)—K'(xa,y»)]  for  both 
subsonic  and  sonic  flow  are  given  in  appendix  C  together 
with  some  remarks  on  evaluation  of  the  kernel  function. 

TREATMINT  OF  TBS  SONIC  CASS 

Because  of  its  special  nature,  the  borderline  case,  Af=l, 
between  subsonic  and  supersonic  flow  deserves  and  requires 
separate  treatment. 

As  Af-»1,  the  e.xpression  for  the  kernel  function  given  in 
equation  (20)  becomes  indeterminate.  It  is  possible,  how¬ 
ever,  to  obtain  conditional  limiting  values  for  the  kernel  by 
considering  the  integral  F,  equation  (4),  and  breaking  it  into 
two  integrals,  Fi  and  Ft,  as  was  done  for  the  general  case. 

With  regard  to  Fi,  its  limiting  value  and  the  value  of  its 
derivatives  with  respect  to  z  at  z=0  can  be  shown  to  be  zero 
as  Af-^1.  From  the  form  of  Fi  given  by  equation  (14), 


lira  F,  =  i|m|ii:.(^Vy?+?)-7[/o(f 

,  ,..oo.r(sv5?+?)'i 

»  •  - - ^dr 

Jo  Vl"i"r* 


Jo  yi  +  T* 
But  since  (see  ref.  22,  p.  172) 


-r^^dr=-Kdii:) 

Jo  Vl  +  r* 
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and  (see  ref.  22,  p.  ;i;i2) 

;  f  m)-Lm  (^5) 

Jj  ^l  +  T-  - 

it  may  be  euucluded  from  e(|uation  that 

lim /’i  =  lim(^^'^=0  {;{6) 

\t-\  .0-1  \  Or  / 

'Pile  total  eoiUrilmlion  to  #f(Xo,i/ol  si  .l/=l,  tlicrefore,  arises 
from  tl\e  limit  of  t',,  equation  (ti),  us  The  limiting; 

form  of  Fi  may  he  written  in  terms  of  nondimensionol  eo- 
ordinates  as 

limf]=:lim|  - — ^  ~= — dX  (.17) 


in  approaching  the  limit  .IP^l  (from  the  subsonic  side)  in 
equation  (37),  it  is  convenient  to  replace  d/by 

d/=l-€ 

where  t  is  infinitesimally  small  so  that 

/S*=(l-.l/)  (l+d/)=«(2-«)«2« 

With  this  appro.vimation,  equation  (.37)  may  be  written  as 

”Jo  - - X - (forz,>0)  (38) 

From  physical  considerations,  the  right  side  of  equation  (38) 
is  to  be  considered  zero  for  Zo^O.  This  is  in  keeping  with 
results  that  would  bo  obtained  if  the  limit  under  consider¬ 
ation  were  sought  from  theory  of  supersonic  flow,  d/>l. 

The  integral  in  equation  (38)  cannot  be  completely 
expressed  in  terms  of  known  functions.  Furthermore,  since 
it  is  singular  at  its  lower  limit,  further  treatment  is  required 
to  reduce  it  to  a  form  such  that  its  derivatives  with  respect 
to  z  can  be  numerically  evaluated.  For  this  purpose  the 
integral  may  be  written  as  two  integrals,  namely 

(39) 


r’  ® 

f  v»**+^ 


The  limits  of  integration  in  equation  (40)  are  so  chosen  that 
the  integral  in  this  equation  can  be  reducctl  to  a  known  form 
by  making  the  substitution 


-r'-J' 

Jo  ^ 


J-»  ,-.UV»i‘+e)r 

,  -  dr  (42) 
0  .y  1  r" 


Equation  (42)  may  be  written  in  tiu-ms  of  the  integrals 
involved  in  Ft  (see  eqs.  (34)  and  (35)),  namely, 

(43) 

Diirercntiating  this  result  twice  with  respect  to  ::  and  then 
setting  c=t)  gives 

Diiferentiating  equation  (41)  twice  with  respect  to  z  and 
setting  z=0  gives 

el  1  ,  i  f"  « 1  (451 

After  perfonning  an  integration  by  parts  and  collecting 
terms,  equation  (45)  may  be  written  as 

<«> 

Equations  (44)  and  (40)  are  combined  to  give  ’ 

Then,  in  accordance  with  equation  (4),  there  is  obtained 
for  /iC(zo,yo)jr.i; 

For  Zo>0, 

iiC(xo,yo)w.i=|r«"’'“"  ■{~/;|^(  ^•(*l2/oi)-2i|^ 


and,  for  ZogO, 


^(*o.!/o).w.i=0 


The  integral  appearing  in  e(|uation  (47a)  is  finite  and 
proper  anti  can  be  evaluatetl  l)y  tiumericul  procedures. 

TRBATMENT  OP  THE  STBAUV  AND  INCOMPRESSIBLE  CASES 

It  is  of  interest  to  consider  the  form  of  the  kernel  function 
given  in  equation  (20)  for  some  particular  values  of  .1/  and  k. 
In  the  following  sections  a  discussion  is  given  for  the  steady 
I'ose  (it=0)  and  the  incompressible  case  (.\/=0).  The  two- 
dimensional  case  is  handled  in  appundi.x  B. 


39461S— 5« - s 
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Beduction  of  the  kernel  for  ‘i  e  cnee  of  steady  flow.— In 
order  to  obtain  tlie  reduet'oe  of  (lu>  kernel  for  tlie  ease  of 
sleatly  How,  eonsider  ilie  expanded  form  ^iven  by  eipiutiou 
{261.  -Vs  k-*0.  there  lesidta  the  followiiu:  expression 

+  -  (48) 

'  \yo  'Ja\jrn+^yn/ 

will'  h  represents  the  downwa.sh  of  a  pressure  <loul>l(‘t  for 
steady  How.  Tliis  result  serves  us  a  partial  idiiK-k  us  to  the 
eorreetness  of  the  expression  for  K(/a..Vo)  ijiven  bv  e»|Utttu»n 
(20). 

Hy  replttiiii}:  j/o  in  <>(|Uttliou  (4.S)  by  ij  —  ti  and  inie^rntin;; 
from  —1  to  I  with  respeet  to  n.  there  is  obtained 


»■/-  ..^J  ifi'o+xJto'+d'fy— I)'  xoH-\/o*+)#'(y+i)'l 

K(xo.Wd,=  -yL— - ^7(7+"!)“  ■  J 


where  the  symbol  ^  indicates  that  a  principal  value  or 

Unite  part  of  the  improper  inteirral  must  be  taken.  (.See. 
for  example,  ref.  24  for  a  discussion  of  Unite  parts  of  such 
integrals.)  This  result  corresponds  to  the  downwash  pro¬ 
duced  by  a  simple  horseshoe  vortex  two  units  wide.  An 
equivalent  expression  for  incompressible  How  is  given,  for 
example,  in  reference  25,  where  in  contrast  to  the  present 
notation,  Xq  has  been  chosen  as  positive  forward. 

Baduetion  of  the  kernel  for  M=s0. — Itt  order  to  etTe<a  the 
reduction  of  the  kernel  for  the  incompressible  case,  the 
expressions  for  Fi,  equation  (15),  and  F2,  equation  (18),  will 
be  examined  for  the  limit  A/-»0; 

From  equation  (15) 

Ira  {-«’<(*lyoi)-y  (A«:ly,|)-/,(iHy.l)l+i}  (50) 

and  from  equation  (18) 

ijt  . . . 


f. 

H  a2’"lyolJo 


i*»l  siuh  •  f/tf— 


yoV»)’+yo* 

Integrating  by  parts  yields 


lim  ^^’=■-,--+•4  - 

"MuS  wtS-fo'+s/o' 


Combining  the  results  from  F,  and  F2  gives  for  the  kernel 
function 

yHxo^+yo-  yo 


•  I  \  X'-(- »/„•  I “*i/X  r  (."):)) 

y«  jii  j 

11 V  .setting  x„  =  l)  in  e(piation  i-'):!).  a  form  is  oblainetl  wliidi 
<-ttn  be  shown  to  agree  with  results  ileriveil  by  KOssner  for 
the  etts<'  .l/=(),  x,i  =  l)  (ref.  (2l>)). 

A  SBUES  EXPANSION  WITH  BESPECT  TO  "A” 

All  approximation  for  tlie  function 

IJffro.j/,,)— A''(r„,i/o)| 

for  small  values  of  k  can  be  olitniiied  by  making  use  of  the 
s«‘ric>s  expansions  for  A',  (eq.  (2.i))  and  for  (/,—£,)  (eq.  (24)) 
and  expaiidiug  all  other  tenr.s  of  A(r„.i/o)  (e<|.  (20))  into  a 
|>ower  series  in  tenrs  of  k.  .\fter  perfomiing  thesi*  expan¬ 
sions  and  eolleeling  tenns  with  respect  to  powers  of  k.  there 
is  obtained  for  .\/<l 

\  yoSxn--rd-V  \Xn‘+^y<i 
-T"  \  I'vd* 


^  12APX.+ - - -J+ 


j~[^(12-U/3*-20W+l5/S»-12fl^)yo'-:i2Af*r,‘+ 
4(.W^-f  tU/»- 1  )xo»-h  1 2g'(d/«-fZl/»-  l)x.ife» 

I2fl®y.^  log 

.j^^[(i ii)d/--  i)xoSx„-+fl**y- 
4A/’(5+-\/»)xo®^ 

\Xn  -rpyir 

5/S<(3AP— l)vBVV+/f'yo*J  .  .  .  ^  (54) 

For  vului!s  of  the  paraiiieters  that  satisfy  the  following 
inequalitie.s 


^  (xo— .u-vx„-+d'yo'‘)<i 


equation  (54)  yields  results  that  are  correct  to  within  about 
2  percent. 
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Com*spon<lingly  for  .l/=l,  o(|iiatioii  (47)  fan  be  (‘xpnixIvU 


[o  obtain 

K(a,1/o)v. 

fl 

•  iktri 

« 

3^0  JQ^  3ye*  l/n'e  I 

L  2  4j-o  4xo’ 

96 

9yo*-6Ty«’-()V«^ 

-  -yo 

/ 

“"I 

i\c 

96 

Xo‘  xo’  3j«.Vo- ,  V  Jfo‘  yo'<  M  1(56) 
^20iy"*3‘^  2  "^xo  f2x,''^20xo‘JJ 

For  viilufs  of  till*  |)uranu‘li>rs  tliut  satisfy  tin*  following 
tnt>(|ualily: 

i*<|uaiioii  (ijti)  yit'lils  ri'siilts  that  art*  coiTi'ct  to  within  about 
2  porrrnt. 

'■’ONCLUOING  REMARKS 

Tltc  main  purposi*  of  this  report  was  to  prt*8t*nt  the  kernel 
function  of  the  intetiral  equation  relating;  the  ilownwash  to 
the  lift  ilistribution  in  a  form  that  can  be  eomputeil.  This 
purpose  has  been  achieved  by  tJic  presentation  of  tlic  koinel 
in  a  form  given  in  equation  (20).  This  equation  has  been 
converted  to  a  form  more  suitable  for  calculation  by  isolating 
the  singularities  as  sliown  in  equations  (30)  and  (31).  The 
special  case  of  A/=l  is  given  in  ecpiations  (47).  The  forms 
of  the  kernel  function  for  other  limiting  cases,  namely  k=0 
and  .l/=0,  are  given  in  ecpiations  (48)  and  (.'53),  respectively. 

L.\NGLEV  AeR0XAUTIC.\L  L.tnOK.ATORV, 

•V.ATioNAL  Advisory  Com.mittke  for  Aeronautics, 
Langley  Field,  Va.,  Sej^ember  18.  lOSS. 


APPENDIX  A 

DERIVATION  OF  THE  INTEGRAL  EQUATION  THAT  RELATES  THE  DOWNWASH  AND  LIFT  FOR  A  FINITE  WING  BASED  ON 

REFERENCE  21 


III  ki'oping  with  tlie  conci'pts  of  linear  theory,  tlic  wing  is 
fonsiilerctl  as  a  nearly  plane  impenetrable  surface.  r.A‘t  this 
surface  i!^  lie  nearly  in  the  x^-plane.  as  indicated  in  sketch  1 
of  the  hotly  of  tlic  report,  and  let  it  anil  the  r.  y,  2  coordinate 
system  to  which  it  is  referretl  he  assumed  to  move  at  a 
uniform  speetl  I'  in  the  negative  r-dircction.  At  the  same 
time,  let  each  point  of  the  wing  be  assumed  to  undergo 
harmonic  translations  of  small  amplitude  Z^{x,yjl)  at 
circular  frequency  u  and  let  c  represent  velocity  of  sound  in 
the  medium. 

The  problem  for  an  oscillating  wing  consists  in  solving  the 
wave  equation  subject  to  certain  boundary  conditions.  The 
wave  equation  in  rectangular  coordinates  is 

The  independent  variable  ^  in  equation  (Al)  is  regarded 
herein  as  an  acceleration  potential;  as  such  it  is  directly 
proportional  to  a  perturbation  pressure  field  and  is  related 
to  a  velocity  potential  ^  as  follows; 

In  order  to  complete  the  boundary-value  problem  for  the 

wing,  it  is  desirable  to  calculate  the  downwash  w{x,y,z,t)=^ 

associated  with  <(/.  Assuming  this  downwash  to  be  harmonic 
with  regard  to  time  implies  that  both  potentials  0  and  ^  are 
harmonic  with  regard  to  time  and  can  be  written,  therefore,  os 


«^x,y,2,t)=e‘“‘ 


0(j:.y,2)'l 

iTU-y--)  J 


With  these  e.xpressions  for  0  and  equation  (A2)  becomes 
independent  of  time  and  reduces  to  an  ordinary  equation 
with  one  independent  variable,  namely 

(A4) 

This  equation  can  be  integrated  witli  respect  to  x  to  give 


_  e~^  -  — 

0— -jy-J  'l\ 


where  the  lower  limit  of  integration  is  chosen,  for  later 
convenience,  so  os  to  satisfy  the  condition  that  <t>  vanish 


The  bountlary-valuc  prohlcm  for  the  wing  may  now  he 
e.xpressed  mathematically  as  follows;  lender  the  assumption 
of  harmonic  motion  the  differential  equation,  equation  (Al), 
becomes 

In  ortler  to  insure  tangential  flow  at  the  wing  surface,  the 
potential  must  satisfy  the  tiownwash  condition 

where  w  and  Zm  are  amplitudes  of  velocity  and  displacements, 
respectively,  and  arc  assumed  to  he  known  from  the  motion 
of  the  wing.  At  2=0,  the  pressure 

P=—Pi4')i-t  (A8) 

must  he  zero  at  all  points  (jr,y)  off  the  wing.  At  all  points 
on  the  wing  ^  is  allowed  to  bo  discontinuous  and  the  value 
of  p  at  a  given  point  is  ilcterminctl  by  the  magnitude  of  the 
discon tinuit}’  in  0  at  the  point.  In  the  neighborhood  of  the 
trailing  eilge,  p  must  go  to  zero,  corresponding  to  the  Kutta 
condition. 

One  other  condition,  that  0  vanish  far  ahead  of  the  wing, 
is  inherently  satisfied  by  the  relation  between  0  and  ^  given 
in  equation  (A5). 

The  potential  0o  at  point  (/.  y,  z)  due  to  a  harmonically 
pidsating  doublet  locateil  in  the  j-j/-planc  at  ({,  a,  01  that 
satisfies  equation  (A6)  is 


where 


5 

dz  It' 


/f'=v'(i-«)’+^(y-n)’+d‘r' 


and  the  factor  ..I  is  a  strength  and  dimensionality  factor  that 
makes  possible  different  uses  and  interpretations  of  the 
potential  00-  If  00  is  considered  as  an  acceleration  potential 
and  substituted  into  equation  (A5),  there  is  obtained  a 
corrcspcniling  velocity  potential  0a  which  may  bo  written  os 


as  X—*—  « . 


0o=A  e 


ON  THE  KERNEL  FUNCTION  FOR  FINITE  WINGS  IN  SUBSONIC  FLOW 


11 


where 


{y—ny+0‘z‘ 


The  downwash 


dz 


associated  with  may  be  written  os 


J..  ,x>+? 


(lx 


(All) 


where  x<,—x—i,  w=«/V'/3*,  and  r=(Jv(s/— With  the 
use  of  this  equation  and  the  concept  of  solving  linear 
boundary-value  problems  by  means  of  superposition  of 
elementary  solutions  to  the  governing  differential  equation, 
the  boundary-value  problem  under  discussion  can  be  written 
as  an  integral  equation,  namely 


rr  -lita  a*  r*o  i>c(x-.wv»»+n) 
w {x,y)=lim^ /ij J  L  e  ^  (/{  dij  —  dX 

^  (A12) 

where  S  represents  the  surface  of  the  wing  and  £({,  i;) 
represents  an  unknown  lift  distribution  or  doublet  strength 
on  S.  Equation  (Al2)  may  be  seen  to  correspond  essentially 
to  equations  (1)  and  (2). 

If  the  distribution  function  L{(,  ri)  in  equation  (A12)  is 
determined  in  accordance  with  the  boundary  conditions 
discussed  in  the  preceding  paragraph,  equation  (A12)  can 
be  considered  as  a  complete  solution  to  the  boundary-value 
problem  for  an  oscillating  finite  wing  in  compressible  flow. 
It  is  also  to  be  noted  that  equation  (A12)  can  be  considered 
to  represent  a  solution  to  the  so-called  “indirect"  problem, 
that  is,  that  of  finding  the  downwash  distribution  associated 
with  a  given  lift  distribution. 


APPENDIX  B 


REDUCTION  OF  THE  KERNEL  FUNCTION  FOR  THREE-DIMENSIONAL  FLOW  TO  THAT  FOR  TWO-DIMENSIONAL  FLOW 


The  purpose  of  this  appendix  is  to  show  that  integration  of 
the  kernel  function  K(xo,ya)  from  —  «  to  -1-  <»  with  respect 
to  if=y—ya  leads  to  a  known  result  for  two-dimensional  flow. 
The  kernel  is  first  modified  to  a  form  that,  for  the  present 
case,  is  easier  to  handle.  Then,  after  performing  an  integra¬ 
tion  by  parts  on  the  modified  kernel,  the  form  of  the  kernel 
for  the  two-dimensional  case  is  given  (eq.  (BlS)).  In  addi¬ 
tion,  the  special  cases  of  d/=l  (eq.  (B23))  and  M—0  (eq. 
(B30)}  ore  also  shown. 

The  integration  under  consideration  with  respect  to  q  is 
equivalent  to  an  integration  with  respect  to  Vo,  namely 

ij  K  (X(,,y—v)dfi=lj  A(ro,yo)dyo  (Bl) 

It  is  remarked  in  advance  that  since  z  has  been  made  zero 
in  the  expression  for  K{xo,yo),  equation  (20),  it  is  necessary 
to  employ  the  concept  of  “finite  parts  of  infinite  integrals" 
when  integrating  this  function  across  the  singularities  at 
ya=0.  Use  of  this  concept  gives  the  same  results  that 
could  be  obtained  by  the  more  arduous  task  of  performing 
the  integrations  before  setting  z  equal  to  zero. 

ModifleatioB  of  the  kernel. — In  order  to  effect  the  desired 
modification  of  the  expression  for  K(x„,yo)  given  by  equation 
(20),  consider  the  first  integral  of  the  e.Ypression,  namely 

-P  <“*'■'"1'  dr  (B2) 

This  integral  can  be  written  as 

lim-P  r"VH^  dr-f-P  f  -^1+? 

•-ho  Jo  Jum  _ 

(B3) 

but  according  to  page  .331  of  reference  22 

rt(S4-iijyol)]  (B4) 


whore  lit  is  the  unmodified  Struve  function  of  first  order  and 
Fi  is  the  Bessel  function  of  the  second  kind  of  first  order.  In 
the  limit  as  these  expressions  have  the  following  values: 
For  the  first .  ssion  in  the  bracket  (see  ref.  22,  p.  329) 

lim  //i(«+(*ll/ol)=f/i(ffc|yo|)  =  -I,(fc|yo|)  (B5) 

and  for  the  second  expression  (see  ref.  22,  pp.  77  and  78) 

lim  F,(64-f*|yo|)  =  -^t///”(i*|yo|)+i«/i(ii|yol) 

=y'ii:.(Jfc(y.|)-/,(ifc(y.()  (Bfi) 


where  //i  denotes  the  Hankel  function  of  the  first  kind  of 
first  order.  With  the  use  of  equations  (B3)  to  (Bfi),  expres¬ 
sion  (B2)  ran  be  wTitten  as 


Vl  +  T*  e- “I'of  dT=P  f’  >1  +  e-'*!'"!'  dr+ 
Jutt 


l^-^i:,(fc|yo|)-l-^I/.(iHyo|)-X.(i|y.|)lj.(B7) 

Substituting  this  result  into  equation  (20)  of  the  text  gives 
the  modified  form  of  K{Xf„y„)  sought,  namely 


*(*o,yo)= 


e 


f  I 

^Myy 


e~  «  - 


MyoVV+dV 

(B8) 


Integration  of  modified  kernel. — Since  the  expression  for 
S{x^,ya)  is  symmetrical  with  respect  to  yo,  that  is,  fiTfzt,,— yo)  = 
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/f(^oi  +  !/o)i  ihc  integration  under  consideration  can  he 
expressed  as 

ij  K{xt,\ya\)dy„=-2lj^  K(Xf,,yt)dyo  (Bft) 

where,  on  the  right,  the  absolute-value  signs  on  yo  can  bo 
dropped. 


After  performing  an  integration  by  parts  by  letting 


I  ya  I  yo 


(BlO) 


and 


I.Vfty. 


artw 


9‘ 

**  f*'  J(X-MVS+W) 


Xr/fl  ^  d\  (Bl  1) 


or 

dtt= 


i 


/S*yoXo  ,  ilcMxdyo  , 

11^2  2  I 


tilyo 


(zo’+d^oT’  a;o*+^yo*  V*o*+/3^o’ 


there  is  obtained  for  uv 

i-sL? '  +H' 


(B12) 


'I 


/  1  1  ^  jB^.-VVit'-MWlx 

Mr 


J-W/fl 


1  +  T*€  '**'<•' f/r+ 


p*^(x-.\/Vx*+4^»r) 


ii-r 

MJo 


"O 

(Bi;i) 


This  e.\pression  vanishes  at  its  upper  limit  yo=  <*>  and  is 
singular  at  its  lower  limit  )/o=0.  However,  by  not  making 
2->0  in  the  derivation  of  K(xo,y<,)  until  after  tliis  stage  is 
reached,  this  singular  value  is  canceled  by  other  terms  that 
have  otherwise  been  dropped.  Thus,  tlie  expression  (Bl3) 
may  be  considered  to  be  zero,  which  is  the  value  of  its  hnitc 
part.  The  integration  under  consideration  is  then  reduced 


to  the  value 


du  which  is 


vdv.=2 


irfr^,^ _ 


ikMxt 

+/SV 


V: 


r  — dr-l- 

'V+/9WJ  J-\i/t^/i+‘r 


,  ,  ^  dX  Vdyo 

0  vX’+/3V  J 
I  Jo  \  \  ^0 

(B14) 

The  terms  of  this  expression  arc  treated  separately  in  the 
next  three  equations: 


First  (see  ref.  22,  p.  180) 
2 


/ikM, 

1  ^ 

\  -®  ^ 

Vxo’+/3*yo’, 

^(n-MVS+W) 


dyo=~  -'*^cosh  'dfl 

- X  W  (^1— (*^^~)] 


(B15) 


second 


.  ■'■‘Xv.VET?  vT+?X 


— 2i*IogJ^ 


(B16) 


an<l  third  (see  ref.  22,  p.  180) 


r~  Pt  ,i5(»-"V^’+^r,')  ,  ...  - 

2i’(  dyo  (  *  -  — ;  — dX=2A*  |  “eO'dX  j 

Jo  Jo  \X’-f-^’yo  J®  J® 


t* , 


.«*  J*  I  €  ^ 


w  P  Jo  Jo  P  Jo 


xt/fc* 


»  •vw+^'yo' 


Vaf//o'«(^|X|)dX 
(B17) 
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Substituting  the  results  in  equations  (B15)  to  (Bl7)  into 
equation  (B14)  gives 

ikSf  I 


I J“  Kix^yijdy, - ^  H, 

H.®(^^)]+^Mog 


1±^+ 

M  ^ 

^  Ho<*>  rfX  j.  (B18) 

This  result  is  a  form  of  the  expression  for  the  kernel  function 
of  Possio’s  integral  equation  rdating  pressure  and  downwash 
for  a  two-dimensional  oscillating  wing  in  subsonic  compres¬ 
sible  flow.  It  checks  the  results  given,  for  example,  in 
reference  27. 

Xsdnetion  of  kernel  for  Ar=l. — ^The  kernel  function  for 
A/=l  may  be  written  as  (see  eq.  (47a)) 

Jf(*ot!to)jr.t=»p  [^-^iWyoD— 

(B19) 

The  second  integral  appearing  in  this  equation  can  be  shown 
to  cancel  several  of  the  terms  so  that  the  kernel  becomes 

(B20) 

so  that  the  kernel  for  the  sonic  case  in  two-dimensional  flow 
may  be  written  as 

aJ  ^  •  — 

fj*  JC(xr),pb)jf_i<fyb=>=— ^-j^2e  *  J*  *^|-dyo— 

««»N«  \ 


dyt)  (B21) 


Integrating  equation  (B21)  by  parts  with  respect  to  poi  re- 
tAinitig  only  finite  parts  of  the  integrated  results,  and  making 
use  of  the  relation 


yields 


it  r-  _!£!:!  \  rt*, 


(B22) 


Finally,  the  kernel  for  the  sonic  case  in  two-dimensional  flow 
may  be  written  os 

if  Ji((*6,yo)*r-.dyo=j-yY<"'**‘ 

(  k  ^  r.  \ 

It  may  be  noted  that  the  integrals  in  this  equation  are  readily 
expressible  in  terms  of  Fresnel  integ  ais 


C(x)=f  cos^<*<ft 


and 


5(i)=“J^  ainj 


m 


Bsdnotion  of  kernel  for  AfsO. — For  M^O  it  is  convenient 
to  modify  the  kernel  function  before  integrating  with  respect 
to  y«.  For  this  purpose  use  is  made  of  the  relation  (see  eq. 
(B7)); 

K,{k\yo\)-^^  [/.(fc|y.|)-i,(ifc|y.|)l 


f  Vw>*+x*«"“^rfx 


.P  r 

w 


(B24) 


and  the  relation 


r  (B25) 

Jo  Ifo  J~*o 


If 
y? 


With  these  relations  the  e.xpi^ion  for  K(xo,yo)tt-«,  equation 
(53),  can  be  written  as 


(B26) 


a/ 
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But 


ifc*  r 

y? 


f  Vyo*+^*« 

J  “-^0 


yo* 


-i*  r 

yo*J-^o- 


.  / 


>Vyo’+x’ 

it  W+yo»  aro,  f" 

J/o’  ^y.Vi?+V  J-« 


(yo'+X>)»'» 


c/X 


(B27) 


therefore, 

/  K(xo,y<i)si-i>=‘~~  p~  J  (B28) 

Integrating  with  respect  to  yo  gives 

/•  2  r"  r"  «“**' 

_/(a^.yoW-.dy.=-j«-“^J^ 

- (y,i+XT* 

(B2fi) 

Integrating  each  integral  in  equation  (B29)  and  retaining 
only  finite  parte  yields 

ij  Jr(z«,yo)i,.ody«— ^ 


2  _ 


(  - 2t  i 

\  ab  Jo 


.  —  dX- 


<■ 


cos  tX 


dX+t 


4rk/  1  .  1 

”Z~  V  2Tixo  2t  * 

t[si  (*x.)+|]J-) 


-^Ci  (tib)+ 


(B30) 


where  Ciftib)  sad  Siftzo)  denote,  respectively,  the  "cosine- 
integral”  and  “sine-integral”  functions  defined  as  follows: 


The  results  in  the  braces  of  equation  (B30)  check  with  results 
given  for  this  case  in  reference  14. 


APPENDS  C 


SOME  BEMABK8  ON  EVALUATION  OF  THE  KEBNEL 
FUNCTION 

Exact  expressions  for  the  kernel  function  £(xo,yo)  are 
given  in  equation  (20)  for  0^Af<l  and  in  equation  (47) 
for  Af=:l.  Corresponding  approximate  forms  are  given  in 
equations  (54)  and  (56). 

Equations  (20)  and  (47)  an  valid  for  any  set  of  values  of 
M,  k,  Xot  and  yo.  To  calculate  the  value  of  the  kernel  from 
these  equations,  it  is  necessary  to  evaluate  numerically 
the  integrals  which  appear.  Vidues  of  the  other  terms  can 
be  obtained  by  making  use  of  existing  tables.  Extensive 
tables  of  the  Bessel  functions  Ki  and  /|  may  be  found  in 
nfennee  28  and  a  table  of  the  Struve  function  Li  with  second 
and  foiuili  differences  for  interpolation  purposes  may  be 
found  in  reference  20.  Sample  values  of  Ae  kernel  are 
given  in  table  I. 

For  certain  ranges  of  values  of  M,  k,  Xo,  and  yo,  as  indicated 
by  equations  (55)  and  (57),  the  kemd  can  be  evaluated 
by  making  use  of  Uie  power  series  expansions  given  by 
equation  (54)  for  0$M<1  and  equation  (56)  for  M=l, 

The  various  expressions  for  Kix9,y»)  become  singular  when 
yt—y-v^iO-  In  order  to  be  able  to  evaluate  the  kernel 
in  such  circumstances,  it  has  been  separated  into  two  parts 
as  shown  in  equation  (30).  One  of  these  is  denoted  by 
K{xt,yi)—R{x^yii  and  is  not  singular;  the  other  is  denoted 


by  A'(2b,yo)  and  contains  all  the  singularities.  Obtaining 
the  value  of  {K—K')  from  the  form  of  the  expression  given 
in  equation  (30),  however,  may  be  troublesome.  This 
particular  value  for  yo=0,  zo>0  can  be  obtained  from  the 
following  limiting  form: 

where  y  denotes  Euler’s  constant  (y=0.577216)  and  Ci  and 
Si  denote  cosine-integral  and  sine-integral  functions,  re¬ 
spectively.  (These  functions  are  tabulated  in  reference  30). 
For  M—\,  this  expression  reduces  to 

/T  «-“*0  file  ik 

lim  (A(r,,y»)-K'(zo,yq)J=--^.{  ^  e »  -|- 

iEC?  t 

f[|-»r-lo,(^)+Ci(^)+iSi(^)-=]}  <C2) 
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The  kernel  function  is  not  singular  for  xo<^0.  For  yo=(i 
and  Xi)<0  it  may  bo  written  for  M<C  I  as 

The  expression  tor  K—K'  for  xoKO,  y«=0  may  also  be  useful. 
It  is 

-^',0  '  Urn 

<«(S)+?]}  ®‘) 

For  K(—x<t,yti=K(—x«,y»)—K\—Xt,y(t)»0. 

Some  results  of  evaluating  the  kernel  and  its  nonsiagular 
part  are  given  as  examples  in  table  I.  (In  order  to  obtain 
these  results  the  required  integrations  were  performed 
numerically  by  manual  computing  methods.) 

APPEimiXD 

ALTEBNATE  POEM  OP  EQUATION  (20) 

Subsequent  to  the  derivation  of  equation  (20)  as  given  in 
the  text,  it  was  found  that  the  two  integrals  involved  in  this 
equation  can  be  combined  in  a  manner  that  leads  to  a  more 
concise  and,  for  many  purposes,  a  more  convenient  form  of 
expression  for  the  kernel  function.  The  purpose  of  this 
appendix  is  to  derive  this  alternate  form. 

Consider  first  the  integral 

(D.) 


and  make  the  substitution 


^  - %«|r 


X=ifc|y,|(Jlf^aT?-r) 
This  substitution  gives  for  Oi 

wi+f 


r-ma  /  Mt  \ 

Ju»  _  \Vi+T*  / 

1  rj{.^vvswi» 

'  J+ 

-JU 

»  _ f 


Vh^ 


Consider  now  the  integral 


XJjglg  ______ 

Vl  +  'T*  dr 


and  integrate  by  parts  by  letting 
u=.vl+^ 


so  that 


dr 


_  T  dr 

■\  1  +  T‘ 


»=r-| - T  e 

I'lyol 


This  integration  gives  for  Gi. 

■ 

-=!=■  (DB) 

^IVol  kipol  k|pol  Jo  vT+? 

Subtracting  Gi  from  (eqs.  (D4)  and  (D6))  gives 

,•  litfuoj  ,  lUHii 

~  t _ 


klPol  /Sklpol 


Mk*yo*' 


1  s  (n-MVn'+S’W) 


SiF5o?‘ 


-fd==  e-^‘‘'»«'dr+ 
%olJo  vTf? 


-a^  (n-AfV5+w) 


^/TT? 


,-*mrdr 


i  iMk\yo\-hfi  1  Sfo-srySHwg) 

Af/u(kpo/  * 

Substituting  this  result  into  equation  (20)  of  the  text  gives 
for  Kixt,yt) 

if(a^,yo)=^^{-ji-fii:,(k|y,|)-g^(/,(k|^^^ 


*l»«l  (kyo)V(k*o)*+^(kpo)* 


^e>^r*fdry 


mjo  "'j 

The  integral  in  this  equation  is  in  general  more  amenable  to 
numerical  evaluation  than  either  of  the  two  integrals  appear¬ 
ing  in  equation  (20).  Furthermore,  with  this  expression,  it 
is  not  necessary  to  consider  the  incompressible  case  as  a 
spedal  case,  since  no  trouble  arises  in  setting  Af=0.  Simi- 
iaiiy,  for  the  sonic  case  no  trouble  arises  and  this  expression 
gives  for  Xo>0: 

ir(*»,yo)K5i— -jjj^  Ki{k\yo\)- 

:  ,  I  (n»n 

^[W|Po|)-I.(k|p.|)l+,|^-(^  S*  L  U.  J  + 


.  eUWr* 
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TABLE  I.— VALUES  OF  THE  KERNEL  AND  ITS  NONSINGULAR  PART  AT  M-0.7 


Xt 

y* 

k 

K(xt,yt) 

*(*t,9*)-IC'(*i,V«) 

0 

a  125 

0.  1 

-63.  827569-t-  1. 112400t 

0. 144529-.  007824* 

.3 

-63.  801759+  a  290793t 

-.  003441-.  060879* 

.5 

-63.  613049+  5  408465t 

-.  009423-.  192655* 

.7 

-63. 127659+  7.  466762» 

-.  0181 14-.  374807* 

1.0 

-621396691+  IR  445693f 

-.  035609-.  756548* 

1.5 

.  125 

.1 

-128  263912+  19.142811* 

.  141754-.  028841* 

.3 

-114.855158+  55.631898* 

-.  031317-.  056060* 

.5 

-9Z  964383  +  88  829346* 

-.  123447-.  115703* 

.7 

-63  878740+109.  927026* 

-.  283001 -.  133318* 

1.0 

-8  792808+ 128  223964* 

-.  581313+.  022309* 

0 

610 

.  1 

-.019271+  .018639* 

-.  000039-.  006699* 

.3 

+.  007493+  .  020950* 

+.  007793 -.  049064* 

.5 

+.  020861+  .  001546* 

.  036165-.  115145* 

.7 

+.  009570-  .  017888* 

.  095337-.  181254* 

1.0 

-.  018833-  .  006290* 

'  .  306627-.  239670* 

1.5 

6.0 

.  1 

-.  027209+  .  020038* 

-.00906  -.006215* 

.3 

+.  002452+  .  028186* 

-.  005415-.  041401* 

.5 

+.021871+  .013305* 

-.  007432-.  109920* 

.7 

+.  022588-  .  008980* 

-.  026790-.  232786* 

1.0 

-.  004786-  .  022987* 

-.  190134-.  523276* 
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